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In this letter we report measurements of the coupling between Dirac fermion quasiparticles (DFQs) and 
phonons on the (001) surface of the strong topological insulator Bi 2 Se3. While most contemporary inves- 
tigations of this coupling have involved examining the temperature dependence of the DFQ self-energy via 
angle-resolved photoemission spectroscopy (ARPES) measurements, we employ inelastic helium atom scatter- 
ing to explore, for the first time, this coupling from the phonon perspective. Using a Hilbert transform, we are 
able to obtain the imaginary part of the phonon self-energy associated with a dispersive surface phonon branch 
identified in our previous work [1] as having strong interactions with the DFQs. From this imaginary part of the 
self-energy we obtain a branch-specific electron-phonon coupling constant of 0.43, which is stronger than the 
values reported form the ARPES measurements. 

PACS numbers: 63.20.D-, 63.20.K-, 68.49.Bc, 72.10.Di 



Topological insulators (TIs), a recently discovered class of 
materials with insulating bulk but exotic Dirac fermion metal- 
lic surface states [2-6], have become the focus of intense 
research by the condensed matter physics community. The 
strong spin-orbit coupling in TIs leads to a definite helicity 
whereby the spin is locked normal to the wavevector of the 
surface electronic state. A fundamental manifest feature of 
such spin-textured surface states is their robustness against 
spin-independent scattering, an attribute that protects them 
from backscattering and localization [7, 8]. In this sense, the 
surface states should be very stable in TIs. Indeed, ARPES 
and scanning tunneling microscopy (STM) have confirmed 
that the topological surface states are protected even against 
strong perturbations, provided they are spin-independent [9], 
even up to room temperature [10]. Consequently, electron- 
phonon (e-p) interaction should be the dominant scattering 
mechanism for surface Dirac fermions at finite temperatures. 
Hence, the study of e-p coupling in TIs is of exceptional im- 
portance in assessing potential applications such as spintron- 
ics. 

In condensed matter physics the e-p interaction plays a 
dominant role in a myriad of phenomena ranging from electri- 
cal conductivity to superconductivity. As such, it has spawned 
an extensive amount of literature. E-p interaction changes the 
dispersion and the lifetime of both the electronic and phonon 
states in a material. The effect of the e-p coupling on the dis- 
persion and lifetime of the states is contained in the complex 
self-energy E (electron) and II (phonon). The real part E' (IT) 
renormalizes the dispersion, while the imaginary part E" (II") 
accounts for the finite lifetime r of the state arising from the 
interaction. In spectroscopic terms, the linewidth (full width 
at half maximum) 
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transformation, it is sufficient to determine either E' (IT) or 
E" (II"). We should note that when determining the self- 
energy associated with an electronic state one integrates over 
all phonon states, and vice- versa. One can then study the e- 
p coupling and its consequences from the electron or phonon 
perspective. Here we report on our investigation of the e-p 
coupling from measurements of phonon dispersions. To the 
best of our knowledge this presents the first attempt at such an 
approach. 

In this report we employ the dimensionless parameter A to 
quantify the e-p coupling. It is defined by [1 1] 
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where a 2 F(uj), the so called Eliashberg e-p coupling func- 
tion, is a product of an effective e-p coupling a 2 involving 
phonons of energy hu and the phonon density of states F(u). 
Previous attempts to quantify the e-p coupling from phonon 
measurements focused on extracting the e-p matrix element 
[12] rather than A, a topic that we reported on in our previous 
paper [1]. 

Because of recent improvements in energy and momentum 
resolution of ARPES, it is now possible to obtain detailed 
information about the strength of e-p interaction as a func- 
tion of the energy and wavevector k of electronic states. 
Yet, extracting such information on the e-p interaction from 
experimental data is not straightforward; it requires approxi- 
mations that depend on the properties of the system studied. 
By measuring the electron's self-energy in the vicinity of the 
Fermi energy (Ep) as a function of temperature, one can fit 
the ensuing temperature-dependence of lifetime broadening 
r e _ p (e^ k) to the relation [11] 



r e _ p (e,,k;T) = 27rA(e<,k)fc B T, 
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is frequently used. Because the real and imaginary parts of 
the self-energy are related by a Hilbert (or Kramers-Kronig) 



valid at high temperatures, and extract the value of A. Yet, us- 
ing ARPES experiments to provide data capable of identifying 
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contributions of individual phonon modes to the e-p coupling 
remains elusive [13]. 

The strong three-dimensional TI Bi2Se3, which is the sub- 
ject of this report, was found to have a single Dirac cone with 
a Dirac point at the SBZ center [14, 15]. Because of this 
simplicity, it has been extensively studied both experimentally 
and theoretically. ARPES studies have been recently extended 
to cover e-p coupling in Bi2Se3 [16-18]. However, the results 
of the reported studies are contradictory. The ARPES mea- 
surements in Ref. 16 and Ref. 17, using Eq. (3), conclude 
that the e-p coupling in Bi2Se3 is exceptionally weak, with a 
value of the dimensionless e-p coupling parameter A ~ 0.08 
[17]. Yet, Ref. 18, using the same ARPES approach and Eq. 
(3), reports a much stronger coupling with A ~ 0.25. All of 
the aforementioned studies are from the electron perspective, 
where the e-p interaction is integrated over all phonon modes. 
By contrast, effects of the e-p coupling on specific phonon 
modes are not commonly investigated. 

From the phonon perspective, it is preferable to express 
the Eliashberg coupling function in terms of phonon mode 
linewidths and frequencies, and write the e-p coupling param- 
eter A as [19-21] 
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where N(Ep) is the electronic density of states at Ep, V a 
is the primitive cell volume in real space, 7j(q, Wjq) is the 
linewidth of the jth mode at wavevector q in eV, cj jq is its 
frequency and Aj(q, o; jq ) is the mode-specific e-p coupling. 
A is then proportional to the sum of the Aj(q, Ujq) of the in- 
dividual phonon branches averaged over the whole Brillouin 
zone. 

Phonon line broadenings due to e-p coupling are generally 
small, even for superconductors with strong coupling. Such 
small broadenings are very difficult to detect in both neu- 
tron and helium scattering experiments mainly because of un- 
avoidable instrument linewidths of a few me Vs. This limits 
possibilities of extracting e-p broadening contributions from 
the measured phonon linewidth. Moreover, in addition to e- 
p coupling, there are other inherent contributions to phonon 
line broadening, such as phonon-phonon interaction (anhar- 
monicity), phonon-defect scattering and phonon anti-crossing 
with other branches. Yet, if IT can be determined, then it is 
straightforward to obtain II" with the aid of a Hilbert trans- 
form. 

Recently we reported experimental and theoretical results 
of surface phonon dispersions on Bi2Se3 (001) [1]. The most 
prominent features are the absence of the acoustic Rayleigh 
branch, and the appearance of a low energy isotropic convexly 
dispersive surface optical phonon branch, here on denoted by 
/?, with an energy maximum of 7.4 meV. It exhibited a V- 
shaped minimum at approximately 2kp, reflecting a strong 
Kohn anomaly. It is the lowest lying surface phonon branch 



within 2k f- Our theoretical analysis attributed this dispersive 
profile to the renormalization of the surface phonon excita- 
tions by interactions with surface Dirac fermions. The contri- 
bution of the Dirac fermions to this renormalization was de- 
rived in terms of a Coulomb-type perturbation model within 
the random phase approximation (RPA), yielding an effective 
interaction strength [22] for the j3 branch of about 3.4 eV/A 
[1]. In this letter, we show that because the energy renormal- 
ization of the (j branch is mainly due to e-p coupling, we can 
simply use the Hilbert transform to obtain the e-p coupling 
parameter A^ for this specific surface phonon branch f3. 

We carried out inelastic helium atom surface scattering 
(HASS) measurements for in situ cleaved (001) surface of 
Bi2Se3, at different temperatures in the range 80 K - 300 K. 
In Fig. 1 we show the experimental data for 300 K and 100 
K (red and blue dots, respectively) for the (j branch along the 
Y-M direction, as well as theoretical fits obtained with the 
pseudo-charge model and RPA [1]. Although on average the 
energy of the 100 K data points is slightly lower at a given 
wavevector q than those at 300 K, the difference is certainly 
quite small, which justifies carrying out the theoretical calcu- 
lations at T — K. A very good agreement between theory 
and experiment is clearly seen. 
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FIG. 1. Dispersion curve of the isotropic phonon branch f3. The 
red squares represent the experimental data collected at 300 K, while 
the blue circles represent the data collected at 100 K. Vertical error 
bars indicate the energy resolution of our facility. The solid line is 
the result of the computational implementation of the pseudo-charge 
model. The dashed line represents the theoretical RPA calculation of 
the renormalized phonon energy [1]. 



Such measurements were also performed on two different 
samples at fixed incident and scattered angles, while varying 
the temperature. Fig. 2 (a) shows the temperature dependence 
of a phonon mode of the j3 branch at q ~ 0.13 A -1 . The 
phonon energy exhibits a very small increase with tempera- 
ture, which is consistent with the trend shown in Fig. 1. Fig. 
2 (b) shows the temperature dependence of a phonon mode 
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outside the j3 branch at q ~ 0.6 A -1 . Consequently, a linear 
fit to the data has a very small slope, which justifies neglecting 
the temperature dependence, in particular, when we take into 
account an instrument resolution of about 1.5 meV. These re- 
sults demonstrate that contributions from anharmonicity and 
defects can be ignored. It also justifies considering only e-p 
interaction in our theoretical RPA calculations. 
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FIG. 2. The energy of a single phonon event as a function of temper- 
ature for Sample 1 at q ~ 0.13 A -1 (a) and for Sample 2 at q ~ 0.6 
A -1 (b), respectively. The circles represent the experimental data, 
while the solid line stands for a linear fit of the data. The gray band 
represents the energy resolution of our facility. 
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FIG. 3. (a) The real part of the self-energy IT 7 , and (b) its Hilbert 
transform II" ', for the isotropic phonon branch f3. (c) The linewidth 
7/3 = — 211" , and (d) the coupling parameter \p(q) determined from 

7/3- 



Renormalization of phonon frequencies due to the e-p cou- 
pling is described by the Dyson equation 

(hw^) 2 = (tuv { °l) 2 + 2(tuv ( °bRe [II(q, lo^)} , (5) 



where 



and aA°i are the renormalized and bare sur- 



face phonon frequency for the f3 branch, respectively, and 
II(q, cjq^) is the corresponding phonon self-energy. An ex- 
plicit expression for II(q, cj q ^) is given in Ref. 1 [23]. We 
used the best-fit parameters obtained in Ref. 1 to calculate 
iFfaoj^p) for the P branch at different wavevectors along 
the Y-M direction. We obtain the corresponding II" (q, cj q ^) 
using the Hilbert transform 
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We plot the results for U f and II" in panels (a) and (b) of Fig. 
3 respectively. 

The phonon linewidth (FWHM) for the (3 branch, defined 
as [11] 



7 ^(q) = -2n // (q,o; q ^), 
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is plotted in Fig. 3 (c) as a function of the wavevector q along 
the Y-M direction. We note that the shape of 7/3 (q) reflects 
the density of the final states for DFQ scattering. Moreover, 
the asymmetry apparent in the steep drop as q approaches 
2k f ~ 0.2 A -1 stems from the suppression of DFQ scattering 
events connecting states with progressively opposing spins at 
the far sides of the circular Fermi surface, as shown in Fig. 4. 
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FIG. 4. Intra- and Inter-band transitions of DFQs that contribute to 
the renormalization of the prominent surface phonon branch p. q is 
the phonon wave vector. The red arrows indicate the spin helicity on 
the Fermi Surface. 



We obtain the corresponding e-p coupling \p(q), shown in 
Fig. 3 (d), with the aid of the expression 



1 7/3 (q) 

2ttN(E f ) h 2 ul/ 



(8) 



as derived from Eq. (4). Averaging over the function \p(q), 
we obtain an effective e-p coupling constant value for the p 
branch \p = 0.43, which is greater than the reported values 
of 0.25 [18] and 0.08 [17] obtained from the ARPES mea- 
surements. This value of A/3, extracted from our experimen- 
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tal data, presents a lower bound on the actual value of the 
overall e-p coupling constant A. However, additional contri- 
butions from the remaining higher surface phonon branches 
to A are not expected to increase its value significantly above 
A/?. The reason is that the dispersions of each of these higher 
branches appear the same regardless of the presence or ab- 
sence of the DFQs as evidenced from our lattice dynamics 
calculations based on the pseudo-charge model shown in fig- 
ures 2 and 3 in Ref. [1]. Moreover, our result is supported by a 
recent optical infrared study [24], which suggests a strong e-p 
coupling for the 61 cm -1 (or 7.6 meV) optical phonon mode. 
Also, a recent theoretical study based on an isotropic elas- 
tic continuum phonon model, where e-p coupling involves an 
acoustic phonon branch, obtains a value of A ~ 0.42 for thin 
film geometry [25], and A ~ 0.84 for semi-infinite geometry 
[26] of Bi2Se3. We note, however, that our experimental data 
demonstrates the absence of acoustic surface phonon modes 
[1]. 

In summary we use inelastic HASS measurements to deter- 
mine the e-p coupling constant for a specific surface phonon 
branch of Bi2Se3. Direct experimental measurement of the 
e-p contribution to the phonon linewidth is difficult, so in- 
stead we obtain it by extracting the real part of the phonon 
self-energy using an RPA fit to the measured surface phonon 
dispersion curve, and then obtaining the imaginary part, and 
hence the e-p contribution to the phonon linewidth, with the 
aid of a Hilbert transform. This approach is supported by the 
fact that the experimentally measured temperature indepen- 
dence of the dispersion demonstrates that the surface phonon 
frequency renormalization is mainly determined by e-p cou- 
pling. Using this procedure we find an average e-p coupling 
parameter for the specific phonon branch A/3 = 0.43 that is 
greater than the values obtained from the ARPES measure- 
ments. 
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